On freeze-out problem in hydro-kinetic approach to A+A coUisions 
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A new method for evaluating spectra and correlations in 
the hydrodynamic approach is proposed. It is based on an 
analysis of Boltzmann equations (BE) in terms of probabil- 
ities for constituent particles to escape from the interacting 
system. The conditions of applicability of Cooper-Frye freeze- 
out prescription are considered within the method. The re- 
sults are illustrated with a non-relativistic exact solution of 
BE for expanding spherical fireball as well as with approxi- 
mate solutions for ellipsoidally expanding ones. 

Introduction — The hydrodynamic approach to mul- 
tiparticle production in hadron collisions was proposed 
in 1953 by Landau It considers as the initial state 
a very hot and dense thermal gas of particles soon after 
the collision, which then expands hydrodynamically until 
some finite time when the picture of continuous medium 
is destroyed. The latter stage, so called freeze-out, hap- 
pens when the mean free path of particles becomes com- 
parable with the smallest of system's characteristic di- 
mensions: its geometrical size or hydrodynamic length. 
Hydrodynamic models find serious utilization in descrip- 
tion of high-energy nuclear collisions, especially at CERN 
SPS and Brookhaven RHIC. Studying predicted spectra 
of different particle species versus the initial conditions 
(IC) and equation of state (EoS), one could get informa- 
tion on early (partonic) stage of the "Little Bang" , such 
as a possible formation of quark-gluon plasma (QGP) 
or even type of the phase transition between QGP and 
the hadron matter. The problem is, however, whether 
the predicted spectra with given IC and EoS in hydrody- 
namic models are unambiguous. 

The standard and widely used method to get the spec- 
tra is the so-called Cooper-Frye prescription (CFp) dj, 
that treat the system at the decay stage of evolution as a 
locally equilibrated ideal gas at some hypersurface. How- 
ever, this prescription presents some serious problems be- 
cause, usually, the freeze-out hypersurface contains non- 
space-like sectors and, as a result, artificial discontinu- 
ities (shock waves) across those hypersurfaces appear to 
adjust CFp to energy and momentum conservation laws 
Moreover, the results of many studies based on cas- 
cade models contradict the idea of sudden " freeze-out" : 
particles escape from the system during about the whole 
time of its evolution and do not demonstrate the local 
equilibration at late stages. The method of continuous 
emission gave an important step in the description of 
particle freeze-out from 4-volume within hydrodynamic 
approach. However, it is not fully satisfactory because it 
is not based on Boltzmann equations (BE) and, as con- 



sequence, when applied to realistic systems, fails as the 
particle escape probability becomes large. In the present 
paper we propose an approximate method for describing 
spectra within hydro-kinetic approach, which is simpler 
for dense systems than pure microscopic one, and apply 
it to a description of particle production from an expand- 
ing and interacting Boltzmann gas that is initially in local 
equilibrium. 

Boltzmann equation, probability to escape and escaping 
function — The BE for distribution function f{x,p) in 
the case of no external forces has the form: 



dfix,p) 



= F^'™(x,p)-F'°^"(x,p) 



(1) 



The term F^"^'" and F^°^^ are associated with the num- 
ber of particles which respectively came to the point 
(x,p) and leave this point because of collisions. The 
term F^°^'^{x,p) = R{x,p)f{x,p) can easily be expressed 
in terms of the rate of collisions of the particle with 
momentum p, R[x,p) —< avrei > n{x). The term 
pgam j^g^g jjjQj-g complicated integral structure and de- 
pends on the differential cross-section. To develop a 
corresponding method, let us split the distribution func- 
tion at each space-time point into two parts: f{x,p) = 
fintix,p) + fesc{x,p), X = (t,x). The first one, fint{x,p), 
describes the fraction of the system which will continue 
to interact after the time t. The second one, fesc{x,p), 
describes the particles that will never interact after the 
time t. Of course, the latter expresses only probabilities 
and cannot be ignored when one is solving BE: possible 
interactions with the escaping part have also to be taken 
into account, in order to find the escape probabilities. 

The function fesc{x,p) is built up as follows. Let us 
denote by x' = {t' ,x+{p/pQ){t' —t)) the space-time point 
where a particle in x with momentum p would be, if it 
moved freely. Consider, at each vicinity of the phase- 
space point {x,p), the number of particles that have es- 
caped from the interacting system during the time in- 
terval (t',t' + dt'). This additional portion of escaped 
particles can be produced only from the interacting part 
of the system. Also, these particles are only among those 
that came to the phase-space vicinity of the point (x',p) 
just after an interaction during the time dt' , suffering 
the last collision there. Therefore, the additional contri- 
bution to fesc{x,p) from the time interval (t',t' -I- dt') is 
'Pix',p)F3"-'"-{x',p)dt' for t' < t and is zero for t' > t. 
Here V{x,p) is the probability for any given particle at x 
with momentum p not to interact any more, propagating 
freely. So collecting all the contributions starting from 



some initial time , we have 

t 

fescix,p)^ fesc{xo,p)+ J dt'V (x' , p)F^-'" (x^ p) , (2) 
to 

where xq = (to,x + (p/po)(io - i)) and fesc{xa,p) cor- 
responds to the portion of the system, which is aheady 
free at tg . It foUows from (||) that 



1 p'^ d 
V{x,p) p° dx^^ 



Uc{x,p)=FS^^-{x,p). 



(3) 



The escape probabihty P{x,p) can be expressed ex- 
phcitly in terms of the rate of coUisions along the world 
line of the free particle with momentum p as 



/ f 

V{x,p) = exp - / dt'R{x',p) 



(4) 



or it satisfies the differential equation 



1 5 j"'°^^(x,p) 

-P{x,p) ^ R{x,p) ^ —- — . (5) 



f{x,p) 



On the other hand, according to the probability defini- 
tion. 



fesc{x,p) = V{x,p)f{x,p) . 



(6) 



This equation, where fesc is given by eq. (H) and V by 
eq. (^), is one of the integral forms of BE (^. One 
can check directly that / = fesc/V is governed by BE 
and, thus, our definitions of fescix,p) and T'{x,p) are 
consistent with BE. 

For initially finite system with a short-range inter- 
action among particles, actually the system becomes 
free at large enough times tout, so 'P{x,p) 1 and 
fesc{x,p) — > f{x,p) in this limit. Our proposal is to 
utilize, for the description of particle spectra in A-f A col- 
lisions, the escaping function with V{x,p) (and thus 
R{x,p), see eq. (^)) and i^f'**" evaluated just in the lo- 
cal equilibrium (l.eq.) approximation for f(x,p). In this 
case the function / — fesc/'P corresponds to a solution 
of kinetic equation in relaxation time approximation. 



df{x, p) f{x,p)- fl.eq. (x, p) 



dxt^ 



'{x,p) 



(7) 



with the relaxation time t{x,p) being expressed through 
the rate of collisions Ri ,,qXx,p) of a particle with mo- 
mentum p, T — 1/Ri.eq. ■ When t —> oo, the relaxation 
time r — > oo also, indicating a transition to the free 
streaming regime. The eq. (0) can be derived also from 
BE with approximation of r.h.s. by F^"™ = Ri,eq. fi.eq. 
and = Ri.eq. f l^- Note that the macroscopic pa- 

rameters of l.eq. distribution fmiction are governed, in 



general, by non-ideal hydrodynamic equations that are 
derived from eq. (|^) in the standard way (see, e.g., 

The approximation proposed is based on collecting all 
the liberated particles during the system evolution and, 
therefore, takes into account rather non-equilibrium pro- 
cesses, similar to particle emission from the periphery of 
a system. Obviously, it should describe well the spectra 
in a fast transition of the system from an initial equilib- 
rium state to free streaming (e.g., explosion), as well as 
in a fairly smooth transition when l.eq. is nearly main- 
tained till rather low densities. Formally, comparing with 
the exact solution of BE, we neglect in eq. (||) the term 



VPS-'" - Pl.eqRl.eq. fl 



l.eq. 



(8) 



integrated along the world line of a particle from to to 
tout ■ At part of a trajectory crossing sufficiently dense 
space-time region, neither F^"™ nor P differ much from 
the corresponding hydro terms in eq. (H). At the later 
stage of expansion, when the densities are small (nor- 
mally n{x) ^ t^^), the values of F^"'™ terms, both for 
hydro and exact solution, are rather small. So the inte- 
grated contribution of the neglected term, eq. (^ , com- 
ing from that part of a trajectory, is also small no matter 
how much the l.eq. distribution function //.eg. deviates 
there from the exact solution /. As for the most delicate 
transition region, one can expect that the neglected term 
(^ in this region could not be important, too. In details, 
the transition from initial l.eq. state to free streaming is 
rather fast at a periphery of the system, so the correspon- 
dent integrated contribution of (||) is small. The contri- 
bution of the neglected term from trajectories crossing a 
central part of the system could be, for l.eq. initial dis- 
tributions, also small if the deviation from l.eq. is slight 
until rather small densities. We shall discuss below some 
examples of such a behavior. 

Particle spectra and correlations — To describe the 
inclusive spectra of particles 

' ^°^° dpt^p2 " Ki^'t'^P^^'P^) ' (9) 

the asymptotic equality fesc{x,p) = f{x,p) can be used, 
replacing the total distribution function / in all irre- 
ducible averages in (||), 



dcr^p^ exp {iqx)f{x,p) , 



(10) 



by fesc ■ Here, p = {pi + p2)/2 , q = pi - p2 and the hy- 
persurface aout just generalizes tout ■ Applying the Gauss 
theorem and recalling that dfi[p'^exp{iqx)] — for parti- 
cles on mass shell, one obtains, using respectively general 
equations (|) and (0) and supposing their analytical con- 
tinuation off mass shell. 



) ^p'' da^fesc{xo,p)e'" 



+ / / °"'d4a;7'(a;,p)i^9""(x,p)e''?^, (11) 

J (Jo 

+ p'^f °d^a;(FS""(x,p)-^^'°^"(a;,p))e*«". (12) 

J (To 

Thus, the use of escaping function as the asymptotic 
interpolation to the solution of BE is equivalent to tak- 
ing, as the source function for the spectra and correla- 
tions, the 4- volume emission function S = -ppaarn 
gether with direct emission fescixo,p) from an initial 
3D hypersurface ctq- The CFp, defined by eq. (|l0| ) 
with substitution (Jout ^ "'/.o. and / —>■ fi.eq., treats 
particle spectra as results of rapid conversion of a l.eq. 
hadron system into a gas of free particles at some hy- 
persurface (Jf o.- Formally, it corresponds to taking the 
cross-section tending to infinity at t < t^f „ (to keep 
system in l.eq. state) and zero beyond to-/ „ • Then 
P(t,x,p) = e{t-tcj ^{y.)) (and so f^^c = at t < i^-j, ^ ), 
and S = PFS^''' = 6(1 - U^ J^))fi,,g, in eq. (0). 

Let us now analyze the particle-spectrum formation 
process based on simple analytical models. It was re- 
cently shown ^ that, if at an initial time t = to = a 
non-relativistic ideal gas has an ellipsoidally symmetric 
Gaussian density distribution and a self-similar velocity 
profile u(a;), then the solution for ideal hydrodynamics 
has the form (C = v — u(a;), v = p/m, V = Hi^i 

/(^^x,v).-(^).exp(-^-5:^) (13) 

where 

T / T/ \ ■^^'^ V ■ 

X,X,^-, T = To{^] , = ^x, . (14) 

In the spherically symmetric case, when all Xi arc 
equal, the l.eq. distribution function ( p^ ) with ( p^ ) is 
an exact solution of BE |Q. One can easily check that 
the momentum spectrum as well as the interferometry ra- 
dius, computed at any time t (at any isotherm) or at any 
other hypersurface enclosing the system, are identical to 
those calculated at the initial time t ~ 0. The reason is 
that, in the spherically symmetric case, the l.eq. function 
( |l3| ) makes the l.h.s. of BE (|^) zero, as well as the r.h.s., 
rendering the evolution of the interacting gas similar to 
a free streaming. Then the volume integral in eq. ([l^ ) 
vanishes and in this case CFp formally gives the correct 
spectrum but its physical meaning is completely different 
from the naively expected one. There is no clearly de- 
fined unique freeze-out hypersurface: neither the escape 
probability V nor the emission function S reveal a sharp 
behavior in {t, x) plane, but are rather smooth functions. 
This is demonstrated in Fig. 1 with the average (in v) 




FIG. 1. The space-time {t,r) dependence of the emission 
function averaged over momenta for an expanding spherically 
symmetric fireball containing 400 particles with mass 1 GeV 
and with cross section a = 40 mb, initially at rest and local- 
ized with Gaussian radius parameter i? = 7 fm and tempera- 
ture To = 0.130 GeV. 

emission function < 5" >=< pi^s''™ >=< PRf where 
V is calculated according to eq. (^). 

Therefore, we can conclude that the difference between 
final spectra (and interferometry radii) at aout and what 
could be found at cto is due to dissipative effects (devia- 
tions from l.eq.), which make the integral over 4-volume 
in eq. ( p^ ) non vanishing. The contribution of dissi- 
pative effects can be essential even if the evolution of 
the system is governed with good accuracy by ideal fluid 
hydrodynamics, that take place for fairly high densities 
or/and cross-sections, and so f{x,p) « fi.eq.{x,p). In 
this case, p^df,f{x,p) « p^'df,fi,eq.{x,p) ~ Kfi,eq.{x,p), 
where k does not depend on the density and the cross 
section but tied, in particular, with symmetry properties 
of the hydrodynamic expansion. For the above discussed 
exact spherically symmetric solution of BE, k — 0, but 
in the general case k ^ and, as a result, for high den- 
sities particle rescattering can lead to a serious changing 
of momentum spectra as compared to the initial ones, 
depending on the initial conditions. 

Typically, however, hadron system is not at local equi- 
librium during fairly long later stage of evolution and / 
could largely deviate from fi ,,q. ■ The proposed method 
of escape probabilities with l.eq. approximation for 
V{x,p) and F^"-"'' allows to use hydrodynamic approach 
for calculations of spectra even in the situation when fi- 
nite inhomogeneous systems go through all stages; from 
local equilibrium to free streaming. Let us illustrate this, 
based on the exact solution ( p4[ ) of ideal hydrodynamics 
for an initially compressed ellipsoid in the longitudinal 
direction, containing the same gas as in the previous ex- 
ample. The results are presented in Fig. 2, where single- 
particle spectra at i = 8 fm for particles that "escaped" 
(became free according to eq.(^) up to this moment are 
shown versus those corresponding to frozen-out ones ac- 
cording to CFp, applied to l.eq. distribution function 
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FIG. 2. The transverse (bottom, at pL =0) and longitudi- 
nal (top, at pt = 0) spectra of particles, escaped until t = 8 
fm/c from an expanding ellipsoidally symmetric fireball of the 
same particles as in Fig.l, initially at rest and localized with 
Gaussian radius parameters Xi = X2 = 7fm, X3 = 0.7 fm 
and temperature To — 0.300 GeV. Dashed lines correspond to 
spectra calculated according to CFp, applied to l.eq. distri- 
bution function at T{t) = 0.063 GeV. 

( p^ ) on the isotherm T{t). One can see that the effec- 
tive temperature of the transverse spectrum calculated 
by using fesc{x,p) is higher than the one given by the 
l.eq. distribution function. This is in agreement with the 
results of partonic cascade algorithm of ref. j|] and also 
with those of continuous emission jj]. The longitudinal 
spectrum shows an opposite tendency. 

The solution of the ideal hydrodynamics ( |l3| ) , (|lj) also 
shows that {Xi /Xi)/{Xj /Xj) 1 with time increas- 
ing, so the velocity field of the expanding system tends 
to a spherically symmetric one. Such a tendency of the 
velocity field is preserved, in the central region, also by 
the solutions of Navier-Stokes equation with the same 
initial conditions, if transport coefficients are calculated 
according to Chapman-Enskog (ChE) method (we used 
a hard-sphere model of interaction). Because of this, the 
deviation from l.eq. in the central part of the system, 
calculated within ChE method, is rather negligible un- 
til the density becomes quite small. Then, there is a 
hope to apply CFp for the soft-momentum spectra, since 
these particles are mainly emitted from there. As for the 
hard-momentum spectra, one expects that such particles 
are mostly radiated from the periphery of the system at 
the early times, because of large hydrodynamic veloci- 
ties and fast transition to free streaming there. So, in a 
rough approximation, one can apply the generalized CFp 
taking into account first, p-dependent hypersurface cr/.o. 
and, second, the deviations on cr/.o. of the distribution 
function from l.eq. due to dissipative effects. 



Conclusion — The proposed method allows to de- 
scribe, in hydrodynamic approach for A + A collisions, 
the evolution of the matter from l.eq. till free streaming. 
The method differs from the continuous emission devel- 
oped in refs. where the central object is the interact- 
ing component fint which is approximated by the l.eq. 
distribution function. However it is possible to show, 
by using an exact solution of BE, that fint is neither 
isotropic nor thermal distribution in its local rest frame, 
even at a moderate value of {V{x,p)). 

It is worth to note that our method, applied to the 
evolution of rarefied hadron gas, overlaps with transport 
models like RQMD. In this aspect, the exact solution 
of BE discussed above is useful as a test of numerical 
cascade algorithms. The advantage of the method of es- 
caped particles reveals when the system is not dilute, 
has more complicated interactions and displays a col- 
lective behaviour. The hydrodynamics is still (or even 
more) suitable for the description of such a system, and 
so our method based on calculations of probabilities for 
constituent particles to escape can be used, considering 
a concrete model for the particle interaction with sur- 
rounding medium. 

The analysis of particle-liberation process demon- 
strates the crucial role of dissipative effects in formation 
of one- and two-particle spectra. These effects should 
be taken into account when experimental data from SPS 
and RHIC are treated in hydrodynamic approach. 
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